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Abstract 

Two concepts of plane waves in anisotropic viscoelastic media are stud¬ 
ied. One of these concepts allows for the use of methods based on the 
theory of complete Bernstein functions. This allows for a deeper study 
of frequency-domain asymptotics of the attenuation function and time- 
domain regularity at the wavefronts. A relation between the direction 
of the energy flux density and the attenuation vector is examined under 
much more general assumptions. 
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complex conjugate of 2 : G C 
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Laplace transform of / 
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vL A+ 
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V • w 
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1 Introduction. 


Real viscoelastic media have a universal property which is usually overlooked 
or not fully exploited. The relaxation modulus of a real viscoelastic medium is 
always a completely monotone (CM) function. We recall here that an infinitely 
differentiable function on K+ is CM if its derivatives satisfy the infinite sequence 
of inequalities 

(_1)"/W(i)>0 fort >0, neNU{0} (1) 

In accordance with its physical meaning it is usually bounded but some theoret¬ 
ical models involve unbounded but locally integrable (LICM) relaxation moduli 
|16] . The LICM property of the relaxation modulus and its consequences for 
creep compliance was studied in detail in m- 

A detailed investigation of wave dispersion and attenuation in viscoelastic 
materials with LICM relaxation modulus in one-dimensional problems can be 
found in [201 0 El 0 |9]. It was shown in [iniis] that the wavenumber vector 
multiplied by the imaginary unit is a complete Bernstein function |19] of the 
variable p := —iw, where w denotes the circular frequency. Integral representa¬ 
tions of the attenuation function (the logarithmic attenuation as a function of 
frequency) and the inverse phase speed were thus derived. The integral repre¬ 
sentation allowed asymptotic analysis of the attenuation and regularity of the 
wave field at the wavefront. 

One-dimensional problems do not exhaust all the practical applications of 
viscoelasticity. Anisotropy is an important property of rocks, bones and bio¬ 
tissues. It is often studied simultaneously with viscous properties of materials. 
It however turns out that the three-dimensional solutions of the anisotropic 
viscoelastic equations of motion cannot be represented by the same classes of 
functions as in one-dimensional problems. It was shown in m that Green’s 
functions in such media have to be expressed in terms of matrix-valued complete 
Bernstein functions (mvCBFs). The last concept was developed in [TO] . 

It is common to begin investigation of wave propagation with a detailed 
study of plane waves, their dispersion, attenuation and polarization. This pa¬ 
per is an attempt to review the plane wave concept in the context of anisotropy 
and viscoelasticity using the methods of m- This paper is dedicated to a few 
selected problems and its objective is to demonstrate the appropriate mathe¬ 
matical tools. 

Combination of anisotropy and viscosity leads to some ambiguities in the 
definition of plane waves. One would expect that plane waves represent inde¬ 
pendent propagation modes. This approach however excludes an application of 
the full power of the theory based on LICM relaxation and the mathematical 
apparatus of complete Bernstein functions. In order to retain the possibility 
of using the complete Bernstein functions it is necessary to work with coupled 
plane wave modes (Section |4]). 

An additional topic investigated here is the relation between the attenuation 
vector (the imaginary part of the complex wave number vector) and the energy 
flux density vector. For this topic we have chosen a larger class of relaxation 
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moduli, namely (matrix-valued) causal positive definite (CPD) functions. While 
the assumption that a relaxation modulus completely monotone is a guess based 
on experimental data, the CPD property can be derived from the fluctuation- 
dissipation theorem using the Wiener-Khintchine theorem. More importantly, 
the theory CPD function is more suitable for the problem at hand. 

The paper contains a brief presentation of matrix-valued complete Bernstein 
functions (Section!^. In the following two sections complete Bernstein functions 
are used to discuss two concepts of plane waves in anisotropic viscoelastic media. 

In Sec ini we define the causal positive definite functions and use their prop¬ 
erties to prove that in a viscoelastic medium with a causal positive definite 
relaxation modulus the angle between the energy flux density and the attenu¬ 
ation vector is acute. This statement is a fortiori true for media with LICM 
relaxation moduli. 


2 Matrix-valued complete Bernstein functions. 

Definition 2.1. A matrix-valued function A : R+ U {0} —>■ Aid is said to 
be a matrix-valued complete Bernstein function (mvCBF) if it has an analytic 
continuation to the cut complex plane C\ ] — oo, 0] which satisfies the inequality 

3z3A(z) > 0 (2) 

and lim^^^o-f -A-(a;) exists and is real. 

By definition 5JA(z) is Hermitian, hence the inequality in ([5]) makes sense. 

Theorem 2.2. If the matrix-valued function A.{x) on]R_|_ is a mvCBF then there 
are two positive semi-definite matrices B and C, a positive Radon measure p 
on R+ satisfying the inequality 



(1 -I- r) ^ A^(dr) < oo. 


(3) 


and a measurable p-almost everywhere bounded and positive semi-definite func¬ 
tion M such that 


A(a;) = B -I- a: C -I- a: / (a: -I- s) ^M(s)/a(ds), a: > 0 

^]0,oo[ 


Proof. If B € and > 0 then for every v G 


(4) 


1 r 


Bv = — v^Bv — viBv 
2i L 


= —v^' (B-B^) V > 0 
2i ^ - 


If A is a mvCBF then S [v^ A{z) v] > 0 for every v G and Sz > 0. 
Furthermore lim 2 ;_j.o-i- A(a;) exists and is real. For real x the matrix A(a:) is 
real, hence it is sufficient to consider v G R'^. It follows that A(a:) v is a 
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CBF for every v S R and therefore there are two non-negative numbers bv,Cv 
and a positive measure /iv satisfying inequality ([3]) such that 

A(x) V = &v + a; Cv + a; / {x + r)~^ /iv(dr) 

•/]0,oo[ 

m- It is clear that B v, where B ;= A(0). As for Cv, 

Cv = lim A(a:) vl 

This proves that the limit on the right exists, hence a symmetric matrix C can 
be defined by polarization 

2v CZ/ AV — Cv Cvr 


It is easy to see that the right-hand side is a linear function of v and w and 
therefore it defines a symmetric positive semi-definite matrix C. 

Let Ao(x) A(a:) — B — a; C, 


By polarization 


Ao(x) 



^v(ds) 


v^Ao(a;)w= / (x-|-s) ^ Av,w(dr) 

J]0,oo[ 

where 2Nv,v,{I) := fiv+v/il) — for every interval I C M+. Av,w(-f) 

can be expressed in terms of Aq w: 


Av,w(/) 


lim — 
e-s-0+ 7T 


Ao(— s -I- k) w 


S — 1£ 


dr 


for every segment I whose ends are continuity points of the measure Av,w 
[19]. The last expression shows existence of a Add-valued measure N such that 
Av,w(d) = N(/) w. N(/) is positive semi-definite because for every v S R+ 

v^N(/)v > 0. 

Now 


N(/) V -b N(/) w ± 2v^ N(/) w = (v ± w)"*" N(/) (v ± w) > 0 

hence 

Iv^ N(/) w| < i [v^ N(/) V + N(/) w] < i (||vf -f ||wf) ^(/) (5) 

where fi{I) := trace[N(/)]. The Radon measure thus defined on IR+ is positive. 
The measures /iv satisfy the inequality (I3|), hence Nv and /i satisfy the same 
inequality. 

By the Radon-Nikodym theorem m there is a /i-almost everywhere bounded 
function M(r) on R+ such that N(dr) = M(r) /i(dr). □ 
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Corollary 2.3. If A{x) is a mvCBF satisfying m 


N(]a,&]) = - lim f 

7Te^0+7]„(,] 


Ao(-s + ie) 


s — le 


ds 


and /i(]a, &]) = trace[N(]a, 6])] for every regular point b > 0 of the measure N 
and every a > 0, where Aq(x) = A(x) — B — x C. 

Equation (jd]) provides an analytic continuation of the function A (a:) to the 
complex plane cut along the negative real semi-axis. 

The function M(s) can be assumed bounded by 1 almost everywhere in the 
sense of measure /r. 

Corollary 2.4. If A[x) satisfies equation o and then 

(i) A(0) = B; 

(ii) C = lima;^.oo A(a:) /or x€lR.+ . 

Proof, (i) follows from equation (g]); 

Concerning (ii), it is sufficient to prove that L(x) := /r(ds) —>■ 

0 for a: —>■ oo. Indeed, for a: > 1, 


(a::)! < f (a;-l-s) V(ds) < f (1 + s) V(ds) 

<j]0,oo[ </]0,oo[ 


< CX) 


on account of (|31). The Lebesgue Dominated Convergence Theorem implies that 
lima;^oo L(a:) = 0, q.e.d. □ 


Lemma 2.5. 


sin(a7T) f°° a;s“ ^ , 
x'~ = - / -ds for X > 0 


n 


a; -|- s 


Proof. Identity (jH]) is equivalent to the identity 

sin(a7T) /■“ s“-i 

X = - / -ds lov X > 0 


71 


X s 


( 6 ) 


(7) 


In order to prove the last identity we use the fact that the Stieltjes transform 
is an iterated Laplace transform and 

poo poo 

/ ds = / z^-^e-^dz = r{a)y^-^ 

Jo Jo 

poo 

r(a) / e“^^ dy = r(a) r(l — a) a:““^ = — - - 

Jo sm(a7T) 

which proves Q. □ 
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(8) 


This suggests an integral definition of the square root 

poo 

B (s I + B)-i s-i/2 ds 

Jo 

provided that B has no non-positive real eigenvalue. 

If / is a CBF then there are two non-negative real numbers a, b and a positive 
Radon measure p on R+ such that 

f{x)=a-\-hx-\- / -p(ds) 

d]o,oo[ s + a; 

Using this formula the function / can be extended to 

f{B)-=al + bB+ [ B(sH-B)^V(ds) (9) 

^]0,c>o[ 

The function x^!"^ is a CBF and equation (|5]) is a particular case of (0). 

In order to prove that the square root B^/^ is a square root in the algebraic 
sense we shall use the following theorem [TB], Theorem 4.1 (5): 

Theorem 2.6. If f-,g and their pointwise product f g are CBFs then {f g){B) = 
/(B) 5 (B) for every matrix B S where /(B) and 5 (B) are defined by 

equation 

The functions f{x) = g{x) = x^!"^ and f g are CBFs, and the dehnition ([ 8 ]) 
has the form of equation ([9]), hence 

B1/2 b1 /2^B (10) 

The two square roots of a matrix B in the algebraic sense are defined as 
solutions Y of the equation 

Y 2 -B = 0 ( 11 ) 

If a matrix B £ has no non-positive real eigenvalue then equation (fTTI) has 
a unique solution with the property that its eigenvalues he in the open right 
half of the complex plane m- This particular solution is called the principal 
square root. 

Lemma 2.7. IfiRB > 0 then the square root Y = B^/^ defined by equation ([5]) 
is the principal square root. 

Proof Equation m implies that Y satisfies equation (ED. The real part of 
the integrand of ([ 8 ]) can be expressed in the form 

i B (B-bsI)”V (B^-bsl)”^ Bt = 

i(Bl'-tsl)"^ [B^(B-bsI) + (Bf-bsl) B] (B-fsI)-i =U'f [B^B-bsKB] U 

where U := (B -|- sl)“^. Hence for s > 0 the real part of the integrand of ([5]) 
is positive and therefore the square root defined by ([ 8 ]) satisfies the inequality 
> 0. This result along with equation ED implies that the right-hand 
side of dl]) is the principal square root. □ 
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Theorem 2.8. If B G and SJB > 0 then the square root defined by 
satisfies the inequality > 0. 

Proof. 

^ B(B + sI)"^ - (Bt+sl)”^ B^ = 

i(sI + Bt)”^ [(sI + B^) B-B1‘(B + sI)] (B + sI)-i = sU1'3BU>0 

where U := (s I+B)~’^. In view of equation ([5]) this implies that 5B^/^ >0. □ 

In view of Definition I2 .1 1 this entails the following important corollary: 

Corollary 2.9. If A{x) is a mvCBF then the principal square root is 

a mvCBF. 

Definition 2.10. A matrix-valued function A : R+ U {0} —>■ Aid is said to he a 
Stieltjes function if it has an analytic continuation to C\ ] — oo, 0] which satisfies 
the inequality 

^z^A{z)<0 ( 12 ) 

and lim 2 :_>o+ -A-(x) exists and is real positive semi-definite. 

If B S is invertible and 3B > 0 then 

3B-1 = — (B'l - B^~i) = —B-i (Bt - B) B^”! < 0 

27Ti ^ ^ 27Ti ^ - 

This proves the following lemma: 

Lemma 2.11. If A{x) is a mvCBF and A{z) is invertible for z ^]—oo,0] then 
A(a;)“^ is a matrix-valued Stieltjes function. 

Theorem 2.12. If the matrix-valued function A{x) is a Stieltjes function then 
there are two positive semi-definite matrices B and C, a positive Radon measure 
fi on R+ satisfying the inequality @ and a measurable p-almost everywhere 
bounded and positive semi-definite matrix-valued function M such that 

A{x)=B+x-^C + [ (a: + s)“^M(s)At(ds) (13) 

^]0,cso[ 

The proof of this theorem, based on the integral representation 
f{x)=a+ — + ( {x + s)~^ fj,{ds) 

^ .^]0,oo[ 

m with a,b > 0 and a positive Radon measure pL satisfying equation ®, is 
analogous to the proof of Theorem 12.21 

Equation m provides an analytic continuation of the function A{x) to the 
complex plane cut along the negative real semi axis. 

Comparison of Theorems 12.21 and 12.121 yields the following corollary 

Corollary 2.13. If A is a mvCBF then x~^ A{x) is a matrix-valued Stieltjes 
function. 
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3 Plane waves in an anisotropic viscoelastic medium. 

We shall consider complex plane wave solutions 

exp {—iu!t + ik • x) a (14) 

with k, a £ C^, of the homogeneous equation of motion 

pu,tt = V ■ [G(t) * Vu^t] t > 0, (15) 

where G(t) Vu,t is Gijrs * dur^sjdt componentwise. 

Substituting (HI in (Ha we obtain an eigenproblem 

[k^ Qijrs{ — iuj)nj Us + poj'^ ar = 0 (16) 

where we have set k = — i/c n with n £ K.^, = 1 and 

Q{p):=pG{p) (17) 

Let Qn denote the rank-2 tensor Qijrs(—'ioj)nj Ug- We shall assume the 
strong ellipticity conditioiHl 

^kirs Ok ttr > 0 for non-zero vectors a, n 

where G°“ = limt_>oo G(t) = limp_j.o Q(p)- Let G°“ denote the matrix n; 

The matrix-valued function Qn(p) is a CBF, hence it is non-decreasing and 
therefore 

Qn(p) > Q“ := G“ > 0 

and the matrix Qn is invertible. Consequently we can define the matrix-valued 
function 

Kn(p) := ppQn(p)~^^^ (18) 

and recast equation (I16p in a simpler form 

[K„(-icc)2-K2l] a = 0 (19) 

A necessary condition for the existence of a non-trivial solution a of (ITOl) is 
the dispersion equation: 

n := det [K„(-iw)2 - l] = 0 (20) 

Equation dT51) is satisfied if either K„(—iaj)a = Ka or Kn(—ia;)a = —Ka. 

This alternative corresponds to two different propagation directions ±n. 

Equation (EOl) is an algebraic equation of degree 3 in A := hence it has 
three roots = Aj(w,n), j = 1,2,3, which are eigenvalues of Kn(—iw)^. Eor 
each root = Xj{oj, n), j £ {1, 2, 3} equation (fTOl) has at least one non-trivial 
solution a £ C^. This solution can be obtained from the expansion of the 
determinant D in terms of the elements of a row and its 2 x 2 minors. Solutions 

^In viscoelastic terminology this condition means that the medium is a viscoelastic solid. 
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a corresponding to numerically different roots Aj (w, n) are linearly independent. 
If a root Aj(a;,n) is double, e.g. A 2 (a;,n) = A 3 (a;,n) for the given values of w 
and n then the number of linearly independent solutions corresponding to the 
double eigenvalue can be either 1 or 2. Summarizing for each solution of the 
dispersion equation the equation of motion (USD has at least one non-trivial 
complex solution (HI with k = ±iAj(a;)^/^n and A satisfying equation (HU). 
The complex conjugate of m is also a solution of (|T5]) . hence m has a non¬ 
trivial plane wave solution 

u(t, x) = 3? [exp (—icjt -I- ik • x) a] (21) 

for each root Aj(uj, n) of £>. 

Theorem 3.1. K„(-) is a mvGBF for every n G 5. 

Proof. Qn(p) is a mvCBF, hence is a mvCBF fCorollarv 12.91) . By 

Corollarv l2.13l p~^ Qn(p)^^^ is a matrix-valued Stieltjes function. By Lemma r2.11l 
its inverse is a CBF. □ 

Since K„(0) = 0, Theorem 12.21 implies that 


'K-nip) =p'Bn+p [ {p + r) ^Mn{r)p.{dr) (22) 

]0,oo[ 

where /r is a positive Radon measure satisfying the inequality (jS]), Mn(r) is 
an S'-valued function bounded /r-almost everywhere by 1 and Bn G S' is posi¬ 
tive semi-definite. The subscript n indicates the parametric dependence on n. 
Corollary 12.41 implies that 


B 


n 


lim [p ^K„(p)] 

p—^oc 


- 1/2 


lim Qn(p) 

p—^■OO 




- 1/2 


where p is considered a real variable, := G\f.ji rikni, is a real symmetric 
matrix. It can be diagonalized and its eigenvalues are non-negative. A non-zero 
eigenvalue 1/cn*^^ of Bn represents the inverse wavefront speed for the corre¬ 
sponding mode and a plane wave front orthogonal to the vector n. 

Equation p2p implies that the matrix-valued function Kn(p) can be ex¬ 
tended to C \ [0, oo[ by analytic continuation and 


Kn(—iw) = —ia;Cn(a;) -b An(a;) 

(23) 

r- ^ ^ T> , f ^Mn(r) 

C^{uj):=B+ 2 

t/]0,oo] ' ^ 

(24) 

\ t \ 2 f Mn('r) ^ 

An(a;):=a; / 

Ajo.oo] ^ 

(25) 
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Both matrix-valued functions are symmetric and positive semi-definite for uj S 
R. Note that Cn(a;) is a non-increasing function of ui with respect to the usual 
order relation for matrices and its tends to the limit B for w — >■ oo. This implies 
that the eigenvalues of the matrix-valued function Cn(w) are non-increasing 
functions of uj bounded from below by the eigenvalues of B. The matrix-valued 
function An(a;) is non-decreasing and it is bounded from above if fj, has a finite 
mass. 

If v(a;) is an eigenvector of Kn(— iw)^ corresponding to the eigenvalue Xj{oJ, n) 
and satisfying the normalization condition vl v = 1 then aj{ijj, n) := JiAj (w, n) = 
v(w)t An(a;)v(w) > 0 and l/cj(w,n) := — 3Aj(a;, n)/a; = v(w)l'Cn(w)v(a;) > 
0. Thus each root Xj of the dispersion equation has the form Aj(a;,n) = 
—iuj/cj{uj, n) -I- aj{uj, n) and the plane wave is given by the formula 


u(t, x) = JR 


^—iujt+iujn-x/cj(uj,n)—aj(uj,n) n-x ^ 


(26) 


where the amplitude a = yv, 7 e C, is a solution of equation (fTOl) for k = Xj. 


4 Plane waves in anisotropic viscoelastic media: 
an alternative approach. 

The function A^ (•, n) is not in general a CBF and therefore the functions l/cj{-, n) 
and aj{-, n) do not have integral representations analogous to (1^ and (051) . We 
shall now try to work around this problem. 

Theorem 4.1. If v is a constant eigenvector of a mvCBF'Y{y) then the cor¬ 
responding eigenvalue is a CBF. 

Proof. Theorem 12.21 implies that 

Y{y) = A-\-yB + y f {y s)“^M(s) ^(ds) 

•j]0,oo[ 

where A, B are positive semi-definite matrices, /i is a positive Radon mea¬ 
sure on R+ satisfying inequality ([3]) and M is positive semi-definite /x-almost 
everywhere. Let v be a unit vector. Then the eigenvalue 

X{y) = Y{y)v = ay by / (1/-f s)“^ m(s) ^(ds) 

-^]0,oo[ 

where a := Av > 0, & := vf Bv > 0 and m(s) := vf M(s) v > 0 /x-almost 
everywhere. Hence X{y) is a CBF. □ 

Corollary 4.2. //Kn(p) has a constant eigenvector a and K(j),n) is the corre¬ 
sponding eigenvalue then k{-, n) is a CBF and k(— iw, n) = —\uj/c{u), n)-|-a(a;, n) 
where l/c{uj, n), a{uj, n) > 0. 

It is clear that l/c(a;,n) is finite, i.e. c(w,n) > 0. 
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Corollary 4.3. Under the hypotheses of Corollary \4-‘A the functions l/c{uj,n) 
and a(a;, n) satisfy equations (1241) and (1251) respectively. 


The function c(a;, n) tends to c°°(n) := l/v’^ Bv as w —oo. 

The hypothesis of Theorem 14.11 is satisfied in for all the plane wave modes 
in isotropic viscoelastic media and for longitudinal waves when the wavefront 
normal n is parallel to a specific direction. In the other cases the eigenvalue 
k(p, n) is not a CBF. It is therefore more convenient to define plane waves in a 
different way: 


u{t, x) = 3? 


icjt—Kn(—icj) n-x 


(27) 


where Kn s defined by equation (11811 and a € is an arbitrary vector. It is easy 
to verify that (EZD satisfies equation m- Such plane-wave solutions appear in 
the analysis of Green’s functions in m- 

We now take advantage of the fact that Kn(p) is a mvCBF and use equa¬ 
tion (1231) which yields the formula 


u(t, x) = 3? 


Cn(w) n-x]—An(t^) n-x ^ 


(28) 


The matrix Cn(a;) plays the role of the inverse phase speed T>(a;,n) while A„ 
plays the role of the attenuation function A{u!, n). 

We now prove that u(t, x) is exponentially attenuated in the direction of its 
propagation. Define the matrix norm 


|A||d := sup 


wt At Aw 


wtw 


(29) 


The norm ||A||d is unitarily invariant, i. e. for every unitary matrix U and V 


||UAV|U=||A|U. 


We can now show that the function 


g-yK„(i^) ^ Cn(w)-y A„(a;) 

is exponentially attenuated if An(a;) > 0 for uj G M+. 

According to Theorem IX.3.11 in [I] every complex d x d matrix A satisfies 
the inequality 

lle^lU < lle^^lU 

Hence 

Let ao(uj) be the smallest eigenvalue of the real symmetric matrix An(a;). If oq 
is the smallest eigenvalue of a real symmetric d x d matrix A with eigenvalues 
Oj and unitary eigenvectors Vj, then 


,-A|| 


= sup 


OUp A 


J2j=i P exp(-2oj-) 

eU I".-1" 


= e 
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Consequently 


||e-^K„(-i<^)||g < g-aoMy 

If An(w) > 0 for w > 0 then equation eiD implies that the plane wave with the 
wavefront normal n is exponentially attenuated with distance. 

For sufficiently small anisotropic dispersion C„(a;) — r(;(a;)l3 = 0[e] and 
attenuation An(a;) = z(a;)0[e], where w{uj) and z{uj) are some real functions, 
the inverse phase function and the attenuation function in (1301) can be approx¬ 
imately disentangled in the form 

exp(icjC„(a;) y - A„(a;) y) = 

d>(?/,a;) exp(iwyC„(w)) exp(-2/An(w)) = 

'i>{y,io) exp(-2/A„(a;)) exp(iw ?/C„(w)) (32) 

where the functions $(?/, w), w) = 1 -(- O [y^] can be expressed in terms of 
exponential functions of nested commutators of C„(w) and An(a;) by the use 
of the Zassenhaus formula M- 

The infinite product is convergent if \y\ ||wC„(w) + An(a;)|| < 0.596705. 

The matrix A„(a;) is symmetric hence it admits a spectral representation with 
eigenvalues ai(a;),n and eigenvectors Wi(a;,n). Logarithmic attenuation rates 
can be made explicit by applying the second line of equation (IMll and expanding 
the amplitude vector in terms of the eigenvectors w^, a = "^^(0;, n), so 

that 

3 

exp(?/[ia;C„(w)-A„(a;)]) a = ^{y, w) exp(iwrC„(a;)) ^ exp(-?/ an.i(w)) Vi Wi{uj, n) 

i=l 

A different expansion can be obtained using the last line of ((32ll . Expand 
the amplitude vector a in terms of three independent eigenvectors wj(w) of 
the matrix Cn(a;) corresponding to the eigenvalues l/cj{uj). This yields an 
expansion of the plane wave in terms of quasi-elastic modes 

3 

exp(?/[ia;C„(w)-A„(w)]) a = '^{y,uj) exp{-y Anjuj)) exp(iu; y/cn,j{uj)) v] w)(a;) 

The quasi-elastic modes appearing under the sum are coupled by the attenuation 
operator exp(—An(w) y) and by the factor 'L(j/,a;). 

The functions $ and are rather difficult to evaluate numerically, but see 
[5] concerning the numerical implementation of the Zassenhaus formula. 

5 Energy flux density of an inhomogeneous plane 
wave. 

Equation 113 implies the energy conservation equation 

^ + W] -h div = 0 
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where the energy flux density is given by the formula 


— ^kl hfc — Gklmn * ^ ^ 


(33) 


and AW/At = akiA\ik,i/At. The energy density functional W for LICM relax¬ 
ation was constructed in and for CPD relaxation in m- 

We shall now calculate the energy flux density of an inhomogeneous plane 
wave u = 5fi [a exp(—icct -|- ik • x)], with a complex amplitude vector a = a^-|-ia^ 
and a complex wave number k = k^ -|- ik^, and then average the flux over the 
period T = 2n/uj. In contrast to the previous sections we do not assume here 
that the vectors k^ and k^ are collinear. This concept is a generalization of the 
plane waves discussed in Section |31 
The time average is denoted by 



Let C = cos {ujt — k^ • x), S' = sin {oot — k^ • x). The particle velocity and the 
stress are given by the expressions 


u = we-‘^'’' [a^C-a^S] 


and 


^kl — Gklmn * '^m,n 



.ran 


(34) 


The identity 



implies that 


CTfei —we 


R,e CimknGklmni, iw) 6 


— icdi-t-ik^-x 


Hence, noting that (C^) = (S^) = 1/2 and {SC) = 0, 

-{^i) = {ukaM)=u:^e-^^"-^ 

^ \ [ ^ S CL d-ui kn Gklrnni^ 1^) “h dm kn Gf^lmni^ 1^) S 


2 


(35) 


and 



-2k^-x 


Im d}^ Gklmn ( kn , I — 1, 2, 3 (56) 
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For comparison with elastic media it is more convenient to express the energy 
flux density in terms of Q{p) = p G(p): 




-2k'x 

2 


Re 


0>kQklmn{ ; 


1 = 1,2,3 (37) 


For an elastic medium Q is constant (the stiffness tensor). 


6 CPD functions and the direction of the energy 
flux density. 


The energy flux density decays exponentially in the direction of the attenuation 
vector. We shall now show that the energy flux density decays exponentially in 
the direction of the flux. To this effect we shall use a much weaker assumption 
about the relaxation modulus than previously. 

Definition 6.1. A locally integrable function / : K.+ <C is said to be causal 
positive definite (CPD) if 

nOCt /*00 /*00 

^ / f{t){4>*(i>){t)dt = Ke / (j){t) / f (s) (j){t — s) ds dt > 0 (38) 

Jo J — oo J 0 


for every square-integrable function (j) with compact support. 

We have used here the notation fit) := /(—t). In [3] CPD functions are 
called functions of positive type. 

This definition is readily generalized to matrix- and operator-valued func¬ 
tions. In particular we are interested in tensor-valued functions. In this context 
tensors are defined as operators on the space S of symmetric d x d complex 
matrices. The space S is endowed with the scalar product (e, f) := J2k i fki- 
The space of symmetric operators on S will be denoted by C. 

Definition 6.2. A locally integrable function tensor-valued function G : IR_|_ —>■ 
C is said to be causal positive definite (CPD) if 


pOO pC 

K / G^jkiiOieij * eki)it)dt = K / 
Jo J-i 


(e(t), G(5) e{t — s)) ds 


/O J—oo L*^o 

for every S-valued square-integrable function e on M with compact support. 


dt >0 
(39) 


We shall need spectral characterizations of CPD functions. 


Theorem 6.3. G is a tensor-valued CPD function on R+ if the function 


Fit) 


Git) t > 0 
0 t<0 


is a tempered distribution and the real part of its Fourier transform N(^) := 
3fiF(f) is positive semi-definite for w S K.. 
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[3], Theorem 16.2.5. 

Since (e, N(^) e) is non-negative for every e G S' (in the distributions sense), 
it is a positive Radon measure. The polarization argument (see the proof of 
Theorem l2.2l) leads to the conclusion that N is a tensor-valued Radon measure. 

Let n = 1, 2 ..., be a sequence of real Schwartz test functions tending 

in the C° norm to the characteristic function xe of a Borel set LI C K. Since 


‘PniO N(0 dC = 2 y [Ho + F(C)+J VniO d^ = 

^ f F(0v5»(C)dC + ^ 




is Hermitian, so is is limit 


N{E)= hm [ <p„(e)N(de) 

n—¥oo J 


for every Borel i? C M. 

Inverting the Fourier transform and noting that F(f) =0 for t < 0 while F(t) 
is real and F(t) = F(t)^ {Fijki{t) = Fkuj{t)), we have 


_ 1 

G{t) = F{t) + F{-t) = —j ^e-‘«‘ [H(dO + H(dOt] 



N(dO 


(40) 


for t > 0. 

Using a version of inequality ([S]) one simplify the above results: 

Theorem 6.4. If the Radon measure n is the trace of N, J := 

1=1 f ^(O^kikiidf) for every (j) G C°, then n > 0 and and there is an es¬ 
sentially bounded C-valued function L : 1R+ C such that N(d^) = L(^)n(d^), 
|L(^)| < 1 for n-almost all ^ G K. 

where |Tp := J^ijki every T G C. We are now ready to state the main 

result of this section: 

Theorem 6.5. Let G be a tensor-valued CPD function. 

The angle between the time-averaged energy flux density ’4' and the attenu¬ 
ation vector is acute. 

Proof. The Fourier transform G(a;) = G(—iw) of a CPD function G satisfies the 
inequality 5RG(w) > 0, hence 


3? 


^k hi Gklmni 1^) dm ^ri 


3? 




I O 


Ok G{-iuj) am kn k\>0 (41) 
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On the other hand 


Gklmni kl kji Q>m — P O/c 


Hence 


5R 


^kklGklmni i^) Ojm ^ri 


3? 


o-k Gki mn ( i^) ^r; 




Gklmni i^) ^ I — 0 


Subtracting equation from equation (j^ yields the inequality 


29 


Gj^i 

mn (-iw) am K 


k^i > 0 


which on account of implies that 

(«') • > 0 


(43) 


q.e.d. □ 

A relation with the LICM relaxation moduli can be established by means of 
the following theorem 

Theorem 6.6. A locally integrable d x d matrix-valued function A.{t) on R_|_ 
satisfying the condition that for each v G the function A(t) v is non¬ 
negative, non-increasing and convex, is CPD. 

cf [3], Theorem 16.3.1. 

Theorem l6.6l can be readily generalized to operator-valued functions such as 
the tensor-valued functions. In particular every tensor-valued LICM function 
satisfies the hypotheses of Theorem 16.61 is CPD. 

Corollary 6.7. If the relaxation modulus is LICM then the angle between the 
time-averaged energy flux density and the attenuation vector is acute. 

The LICM property is defined in terms of an infinite sequence of inequalities 
([T|) for derivatives of all the integer orders. However only the first three of them 
are relevant for the thesis of Corollary 16.71 This is interesting because the CM 
property cannot be verified from raw data while the first three inequalities m 
can. 

Viscoelastic media with CPD relaxation moduli were studied in [12]. 
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